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We study mesons on the lattice with a special focus on excited states. For that purpose we con- 
struct several quark sources with different spatial smearings, including p-waves. These quark 
sources are then combined with the appropiate Dirac structures to form meson interpolators of 
definite spin. We use these operators to construct a cross correlation matrix from which we ex- 
tract ground and excited meson states using the variational method. For the calculations we use 
gauge configurations with Nf = 2 dynamical Clover Wilson fermions provided by the CP-PACS 
collaboration. We show preliminary results for pseudoscalar, scalar, vector and pseudovector 
mesons. 
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1. Introduction 

Ground state spectroscopy in lattice calculations appears to be well understood. Excited state 
spectroscopy, however, is still a challenging task. There are two major problems: First, one has to 
improve the overlap of the interpolating fields with excited states. Excited hadron states include 
both radial and orbital excitations and thus allowing for such excitations should be implemented in 
a lattice calculation to obtain more realistic results. The second problem is finding reliable means to 
disentangle ground and excited states in the hadron spectrum. Our method of choice is a variational 
approach [|TJ, 0J. It not only allows us to extract excited states but also separates physical states from 
ghost contributions in quenched and partially quenched calculations. 

2. Simulation details 

We perform lattice calculations with dynamical gauge configurations with Nf = 2 Clover Wil- 
son Fermions. The configurations were generated by the CP-PACS collaboration, which made them 
publicly available. At the moment there are three lattices available: a 12 3 x 24, a 16 3 x 32, and a 
24 3 x 48. All these lattices have approximately the same spatial volume of about 2.5 fm, but differ 
in their lattice spacing. This enables us to perform a continuum extrapolation. For the preliminary 
results we are presenting here we have used 100 configurations per sea quark mass of the 16 3 x 32 
lattice which has a lattice spacing of 0.1555(17) fm. For further information about the lattices and 
how they were generated, see Ref. [^, f|]. 

For our simulations, we make use of the Chroma software package ^ from the USQCD. 
This package has the advantage that most of the applications needed in lattice QCD are already 
implemented and it can be easily installed on various platforms, including QCDOC. 

3. The Method 

To obtain reliable results for excited states, we follow a procedure which has already proven 
to be very successful in quenched calculations [Q, f7f|. First, we generate several spatially different 
quark sources. In previous studies, we considered a narrow and a wide Gaussian source (in the 
following called n,w, respectively), obtained from gauge covariant smearing, to allow for a node in 
the radial wave function of the quark. Here we also include p-wave sources (p x , p y , p z ) which we 
generate by acting with a covariant derivative on the wide smeared source. This enables us to not 
only explore the possibility of radial excitations, but also orbital ones. In addition to these sources 
we add a local source (L) so that we can examine physical matrix elements. 

We then combine these sources with the appropriate Dirac structures to obtain interpolating 
fields of definite spin. By doing so, we end up with a large number of interpolators for each 
meson channel (e.g., 15 different interpolators in the case of the pseudoscalar meson). From these 
interpolators we construct a matrix of correlators and then apply the variational method which has 
been proposed by Michael and later refined by Liischer and Wolff [Q]. Here, one has to solve 
a generalized eigenvalue problem. This method has some advantages. First, the system has full 
freedom to choose relative contributions of the different interpolators in the diagonalization step. 
Second, this approach can separate ghost contributions in quenched and partially quenched results, 
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as has been shown in Ref. [Qj. In the case where one uses more than one local operator one can 
also extract ratios of couplings to excited states, see Ref. [Q]. 

The masses of the excited states are obtained by looking at the eigenvalues, which to leading 
order behave as 

k®(t) <* e- tM *. (3.1) 
From these eigenvalues we then construct effective mass plots by using 



2 J + 

When we find a plateau in these plots and if the corresponding eigenvectors are steady we conclude 
that the signal of the considered state is disentangled from higher excitations. We then fit the 



eigenvalue to a single exponential according to Eq. (£3JJ). In the following we present the results of 
these fits for different mesons. 



4. Discussion of the Results 

In Fig. [T], we show a collection of effective mass plots for the different meson channels pseu- 
doscalar(PS), scalar(SC), vector(V), and pseudovector(PV). The effective masses are shown only 
for the completely degenerate case 

= = = (4-D 

With our limited statistics, we find good plateaus for most of the ground states. But the results for 
the excited states do not look encouraging. In fact, we are able to obtain a first excited state only in 
the PS and V channels. 

In Fig. ^-||, we plot the results of our fits versus 



1 
K 




(4.2) 



We destinguish three different cases: The case where both valence quark masses are equal to the 
sea quark mass, denoted as SS, the case where only one valence quark has the same mass as the 
sea quark, denoted as SV, and the VV case where both valence quarks have a mass which differs 
from the sea quark mass. In all the figures the vertical line represents 1 / Kent obtained by the CP- 
PACS collaboration using all available configurations. Since our results are still very preliminary 
we perform only linear fits for the chiral extrapolation. 

In the left hand plot of Fig. [2], one can see that we obtain with our limited statistics a result for 
Kern which is comparable to the CP-PACS result for this quantity. 

In all the cases we cannot exclude systematic effects which are due to the limited statistics we 
have so far. It might, for example, be possible, that with larger statistics some of the plateaus in 
Fig. [j] change and a fit interval starting one timeslice later might be more appropriate. Especially 
our excited states with their short plateaus can be strongly affected by such systematic shifts. This 
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Figure 1: Effective masses for pseudoscalar(PS), scalar(SC), vector(V), and pseudovector(PV), shown for 
completely degenerate quark masses. The horizontal lines denote the time intervals of our fits and represent 
the fit results m± G,„, where O m is the statistical error obtained from single elimination jackknife. 




Figure 2: Results for the pseudoscalar meson ground (left plot) and excited state (right plot) versus 1/fc. 
Filled circles denote the SS case, open squares the SV case, and open diamonds the VV case. The vertical 
line is 1/ fQ nf obtained by the CP-PACS collaboration. We also show a naive chiral extrapolation of the SS 
results. Also the experimental value is included. 
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Figure 3: Same as in Fig. 2, but for the scalar meson ground state. 




might also explain why our results for the excited pseudoscalar and vector meson are extrapolating 
to values which are much higher than the experimental ones. 

Apart from this, we also find that the widths of our Gaussian sources depend on the sea quark 
mass. This makes it rather difficult to tune the smearing parameters to obtain widths in physical 
units which are approximately the same. What we also find is that a cross correlation matrix of 
interpolators built from L and n sources gives much better results than using interpolators built 
from n and w sources. This suggests, that our narrow and wide smeared sources are too similar to 
allow to disentagle ground and excited states. 
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Figure 5: Same as in Fig. 2, but for the pseudovector meson ground state. 



5. Summary and outlook 

We have presented preliminary results of our calculation of ground and excited meson states 
from dynamical configurations with Nf = 2 Clover Wilson fermions. For all considered channels 
we find results for the ground states and in the pseudoscalar and vector channels we are even able 
to extract a first excited state. 

We encounter some difficulties to optimize the smearing parameters for our calculation since 
we find that the width of our Gaussian source depends on the sea quark mass. In addition to this, 
we observe that interpolators built from these Gaussian sources are giving worse signals for the 
excited states than interpolators built from one Gaussian and a local source. Possibly, this means 
that the Gaussians used are too much alike. 

To circumvent this problem we are currently running simulations with a different set of sources. 
We use a narrow source with much smaller width and replace our old wide source by the narrow 
source upon which we apply a Laplacian (The Laplacian being a scalar operator does not change 
the quantum numbers of the hadron). 

Our future plans also include performing the same calculations on the other lattices of the 
CP-PACS collaboration to be able to do a continuum extrapolation of our results. Especially the 
couplings to excited states are in that respect of great interest. 
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